


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 


DSpace Repository 


Theses and Dissertations 


1972 


l. Thesis and Dissertation Collection, all items 


Comparison of numerical and hydraulic 
oceanographic prediction models. 


Salas Romer, Leopoldo. 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/16030 


Downloaded from NPS Archive: Calhoun 


ША | DUDLEY 


(I | orar 


http://www.nps.edu/library 





Calhoun is the Naval Postgraduate School's public access digital repository for 

research materials and institutional publications created by the NPS community. 

Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholarty author. 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


COMPARISON OF NUMERICAL AND HYDRAULIC 
OCEANOGRAPHIC PREDICTION MODELS 


Leopoldo Salas Romer 











em California 93940 


— 





ae 
ar 


NAVAL POSTGRADUATE SCI 


lonterey, Galifornia 


e 
TA 


JUL 





COMPARISON Or NUMERICAL AND HYDRAULIC 
OCEANOGRAPHIC PREDICTION MODELS 


by 


Leopoldo Salas Romer 


Thesis Advisor: BBs Thorsecr 
September 1972 


же" 


ow 
A 
4 
j 


Approved for public release; distribution unlimited. 


1152287. 





Comparison of Numerical and Hydraulic 


Oceanoperaphicenrediction Models 


by 


Leopoldo Salas Romer 
Commander, Venezuelan Navy 


Оше: in partial fulfillment. of the 
requirements "Por the- degree oi 


MASTER OF SCIENCE IN OCEANOGRAPHY 


from the 


NAVAL POSTGRADUATE SCHOOL 
ΕΤ Πως... 


w 


\ 





Library 
Naval Postgraduate School 
Monterey, California 93940 


ABSTRACT 


The Hydro-Numerical Prediction Model developed by Hansen 
ορια το san Diego Bay, and the results compared both 
with the hydraulic model and the real data obtained by 
field measurements. This allows one of the few good compar- 
isons between numertcal and -hydreulic models for the 
Fiction of actual conditions. 

The bay was divided into two sections that were run 
Separately in order to obtain the desirable spatial resolu- 
mem. This division required solving the problems of proper 
mene and matching techniques between both portions. The 
solution involved the addition of an appended pseudo-bay 
Ene first section of the model in order to compensate 
The correct tidal prism. The effects of a proposed 
Aena open entrance in the southern part of the bay were 
Enudsed. This resulted in an increase of flushing in the 
Excbern portion of the bay but caused the currents in the 
center of the bay to be small which decreased dispersion 
in the central portion of the bay. In general, both models 
produced similer and reliable results, but there was a 
Esussderable reduction of cost and time with the numerical 


model. 





г | EN 





TABLE OF CONTENTS 


Е ансо -—— 2 сЕ лыша 
ЕЕЕ AA A 
ου o meee A TT, m 
МОО HYDRAULIC MODEL OF SAN DIEGO BAY ---------------- 
A. “DROOREIPTION OF SAN DIECOeBAY----—------------- 

B. PURPOSE OF THE U.S. CORPS 
СЕНЕК Ее МОРЕ ЭШ а-а L 
А ЕСО аа-а 
D. RESULTS OF THE HYDRAULIC MODEL -------------- 
III. THE HANSEN HYDRODYNAMICAL-NUMERICAL MODEL ------- 
П А ТО SOP THE GP ID amo ae nep tato nte 
СО СООРО ОЕША THE MOBEL (1) ----------------- 
er ADDITION OF AN APPENDED PSEUDOCBAYXS---------- 
OE) o ocio ee nee 
MEO MEE oo RR о cose... 
IV. MERITS OF HYDRAULIC AND NUMERICAL MODELS -------- 
ΠΝ... CONCLUSIONS -------- EEE tk 
КО ШИ А NUMERICAL PROGRAMS -------------------~-- 
БОО yS аы вва 
DISTRIBUTION {Π8Π ----------------------------- 
νιν: --------------δ----ε--πεᾱ---------- 


1“ 
al 
T2 
14 
14 


16 
NS 


e 





21. 


LIST OF FIGURES 


Eee IRE CTE oo 


ВЕСЕ ТОШЕ PROTOMBE AND IN THE 
HYDRAULIC MODEL ----------------------------------- 


A. CURRENT AT A-2 U.S. CORPS OF ENGINEER RESULTS 


CUURENT AT B-2 U.S. CORPS OF ENGINEER RESULTS - 


l 


CURRENT AT C-2 U.S. CORPS OF ENGINEER RESULTS 


SCHEME OF THE GRID NET ------------------------ 


B 
A 
Pee CURR NIT De Us) CORES OF ENGINEER RESULTS 
A 
B COMPUTATION GRID ------------------------------ 
A INPUT BOUNDARY AND FREE BOUNDARY 

WITHOUT CONTINUATION -------------------------- 


B. LEVELING OF THE EDGES (LUBRICATED WALL) 
AT BOUNDARIES -----------------------------.-.--- 


C. CONTINUATION OF CURRENT AT THE FREE OPEN 
BOUNDARY WITH THE GRADIENT CORRECTED BY 
Μο πα ο ο eee E nooo 


A. VALUES OF Z ALONG MODEL (1) FOR DIFFERENT 
CONSTANTS IN THE CONTINUATION OF THE FREE 
BOUNDARY AFTER 1 HOUR OF REAL TIME ------------ 

VALUES OF VELOCITY AT M-3 AND N-55 TO 53 IN 

IHE FREE OPEN BOUNDARY FOR DIFFERENT CONSTANTS 

IN THE CONTINUATION AFTER 1 HOUR OF REAL TIME ----- 

APPENDED AREAS IN MODEL (1) ----------------------- 

CIRCULATION IN MODEL (1) WITH THE APPENDED AREAS -- 


П Pee AT NAVY PIER MODEL (1) FOR 
Dito ΑΕ ΕΡΕ Νανα = 


B. HEIGHT AT PIER N. 2 MODEL (1) FOR 
επ! ἍΡ AROS —————————------------ 


A. CURRENT AT A-2 MODEL (1) FOR DIFFERENT 
APPENDED AREAS ------------------------------- = 


t5 


23 
2 
60 
26 
26 


33 
>> 


39 


41 


43 


46 





M. 


In. 


14. 


Lise 


16. 


E. 
18. 


19. 


20. 


20-1. 


el. 


2. 


23. 


B. CURRENT AT C-2 MODEL (1) FOR DIFFERENT 

I A o is 46 
υπ TRANSPORT IN SECTION AA] MODEL (1) ---------- 47 
EOL ME TRANSPORT ON SECTION AA2 MODEL (1) ------ 47 
A. CURRENT AT A=2 MODEL (11) ----------------------- 48 
B. CURRENT AT B-l MODEL (1) ------------.---------- 48 
E IAS РЕКУ ——-------------------- 49 
B. SNET=NOMVNEDSTRANSRORTJNODEL (Л) -—--------------- 49 
A. - HEIGHT AT NAVY PIER MODEL (1) ------------------ 50 
ВО ш απ Ρις М ο ΜΑΡΙ, (1) ------------------ 50 
CURRENT DISTRIBUTION IN MODEL (1) AT TIME = 0 Hr --- 51 
CURRENT DISTRIBUTION IN τ; ΠΠ - 5:5. 
ЕВЕ ВОО ИЕМ MODEL, (1) AT TIME = 9 Hr === 53 
A. CURRENT AT C-2 MODEL (2) FOR DIFFERENT 

ШЕШ ше 4 4 J l 55 
B. CURRENT AT D-2 MODEL (2) FOR DIFFERENT 

а Ио ΠΤ Τ en en 55 
A. HEIGHT AT PIER N. 2 MODEL (2) FOR 

ΙΡ ΙΙ  Ξ------------------- 56 
B. HEIGHT AT SOUTH BAY MODEL (2) FOR 

II O .--------------------- 56 
IN IO merce ____________ 58 
A. CURRENT AT C-2 FOR DIFFERENT MINIMUM DEPTHS ---- 60 
B. HEIGHT AT PIER N. 2 FOR DIFFERENT 

A sr --------------- 60 
νι, 7, ------------ 61 
D. СОЗА UE ρε ο ك‎ 61 
A. HEIGHT AT PIER N. 2 MODEL (2) ------------------ 62 
υπο СОШ BAY MODEL (2) ------------------ 62 





24. 
ED. 
ΓΕ. 


EN. 


BO. 


29. 


BO. 


Bl. 


Eo. 


33. 
34. 
25. 
Bo. 


3T. 


n 


CURRENT DISTRIBUTION IN MODEL (2) AT TIME 3 Hr -- 


ΤΙ η ηππΠΠῚ IMN MODEL (2) AT TIME 9 Hr -- 


CURRENT DISTRIBUTION IN MODELS (1) 
CN DUE UEM" ''.......----------------:-- 


CURRENT DISTRIBUTION IN MODELS (1) 
aad ООШ T E 


CURRENTE DISTRICUCION IN MODELS (1) 
UGC Ole me 2 ος κα ος... 


CURRENT DISTRIBUTION HYDRAULIC MODEL 
ur mM pie Pj O HR ----—------------------ 


CURRENT DISTRIBUTION HYDRAULIC MODEL 
NIMES CUN ЕН ----------------------------- 


Ας EL DR O 3) ου. 
NE ο hui Dp nf ---------------------- 
ο ΠΠ ο ο CAVE MODEL) (3) =============.= 
oa MODEL (3).----------------- 
ρα πμ πες οι MODE, (5) -------------------- 
RULCDNEDERSBRISUTION MODEL (3) —e-2---------2-------- 
BULLENMSDISBISRSUPION MODEL (3) ———. 


CURRENT ES TRIEUTION HYDRAULIC MODEL 
(PROPOSED SECOND ENTRANCE) AT TIME = 12 Hr -------- 


CURRENT DISTRIBUTION HYDRAULIC MODEL 
(PROPOSED SECOND ENTRANCE) AT TIME = 18 Hr -------- 


64 
65 


66 





NY 


km 


TABLE OF SYMBOLS AND ABBREVIATIONS 


Phase wave speed 
Neme” an the model 
Melo in the prototype 


Total difference 


Ur. ,® SES cse ы QU p: 
*rictton Force 


Coriolis parameter 

Proude number 

Becceteration ol gravity 

Gravity in the model 

Gravity in the prototype 

Total depth (h + n) 

Depth of the water at mean sea level 
пев πω CHE model 


Depth in the prototype 
DEEN al U and Verid points 


DEEN aves erid points 

πι πο ο πο ο Cols yan 

Coefficient of horizontal kinematic eddy viscosity 
ΠΣ ΕΤΟΣ | 
Kilometer 

США ЧО Еле length 

Е све length in the model 


Characteristic length in the prototype 





cm 


MEH 
NEH 


ЧЕИС 


ο ο we MU PU 
< e SES 


с 
< 


<< < 


ZU V 


Length of mesh 
Grid co-ordinates (also m,n) 
Meter 


Centimeter 
Near border Erid points 


Hyaroctatick pressure 

Bottom roughness COCO reren L 
(таум гало, Gravity scale 
Veleciry=racio- Velocity scale 
Depth ratio. Vertical scale 
hemor heat Owe HOP’ Zoncal scale 
Time ratio. Time scale 

Deen προ αυτο Discharge scale 
ip pro l JOSE Volume scale 

το μπαμ ο PirpictCion 

Time 

Components of horizontal velocity 
hovel velocity vector 
Characteristic velocity 

Mean u ana v components 

Wind speed component 

Space co-ordinates 

τη BONES 

smoothing parameter 


l- a 





Partial difference 
Specific weight 

Drag coefficient 

sea level anomaly 
Density oiche fluid 
Wind stress 


porbh-amgubar velocity 


Ec pM 
dx oy 


Bottom stress 


Difference 





ACKNOWLEDGEMENTS 


T wish to express my sincere thanks to everyone who 
has assisted in the preparation of this study ио рагі 
mmn TT extend ny gratitude to Dr. T. Laevastu of the 
mni ronmental Prediction Research Facility, to Mr. Sheldon 
Mazanort “01 beet “Nemertieo _aWeather enbealL,.and to Miss 
Egon Raney of the Computer Center for their helpful 
Sugceestions, encouragement, and programming assistance. 
Шау to Dr. Edward B. Thornton, my thesis advisor, my 
ο. νο thanks for his patience, understanding, and 
ο UCbIVe Criticism that contributed greatly to the 


pletion of this project. 


10 





ENTRO UC TON 


КОКО REVIEW 

Hl Study Of Days, estuaries and other Semi-restricted 
sS iS difficult because these masses oi water are subject 
Mco nstant movement with short time irregular periodicities. 
mie important mechanism driving the Circulation these 
Areas generally is the tide. The water is partially re- 
mnewed and mixed during each tidal cycle due to the inter- 
change ol water between such areas and the ocean. The 
shape and depth of the area governs the direction and speed 
e currents generated during this interchange. The 
constant flow of water causes a constant state of non- 
Саат of constituents, temperature and potential 
ον Inside the bays or estuaries. 

The study of such areas by field observations is ex- 
КО леу laborious requiring large expenditures of manpower 
Ena money. Another approach to the solution of the problem 
K r reproduce the physical characteristics of these areas 
EN Dc use of hydraulic models. Hydraulic models have been 
used for many decades and have proven to be a reliable 
method for studying an area although space and operating 
problems make them a relatively expensive tool of 
investigation. 

With the advent of high speed computers, the numerical 


solution of the differential equations used to describe 


rI 





ШО ama de possible. Oceans, seas, gulfs and other 
water masses have been described by the use of this tech- 
waue. It transforms a physical model into a program 

can δε modified to simulate different conditions more 
NES. can be stored for future uses and alterations, is 
ШЕ ет and allows solutions to be obtained at reasonable 
EROS. 

The U. S. Army Corps of Engineers made a study of San 
éso Bay by means of a hydraulic model based on field 
measurements made in early 1967. The results have been 
reported in a Technical Report [10] and published in a 
paper [1] by the U. S. Corps of Engineers. 

The Hydro-Numerical Prediction Model developed by 
Hansen is applied to San Diego Bay, and the results com- 
pared both with the hydraulic model and the data obtained 
po Ice lid measurements. This allows one of the few good 
eomparisons between numerical and hydraulic models for 


mier prediction of actual conditions. 


Eee OBJECTIVES 

The objective of this thesis is to reproduce by means 
of the Hansen hydrodynamic-numerical model the existing 
rL ons in San Diego Bay and compare them with condi- 
vions produced in the hydraulic model operated by the U. 5. 
EsPps of Engineers. This comparison will result in a 
proper evaluation of the merits and limitations of both 


models. 
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The division of the bay into two separate models is 
EMUempted in order to obtain comparable resolution in the 
merical model. This division requires solving problems 
che method of running numerical models in series and of 
Mastering boundary conditions from one section to 


another. 
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IAE MODEL OP SAN DIEGO BAY 


DESCRIPTION OF SAN DIEGO BAY 

ο. σσ Bay 1s located on the southern coast of 
ИШЕ ουπα, U.5.A. The bay is long, curved, has an area 
of approximately 90 sq. km. and a maximum width of 3.1 km. 
"eig. 17]. "Tt-ts-cenmecebted-wdéth the-ocean by the narrow 
mica Channel at its northern end. 

ОО low oí fresh water can be considered negligible 
τ Curing Focal heavy rainfall. The circulation of the 
entire bay is driven by the tide at its unique entrance. 

a dal je lreulation can be perturbed by local wina 
eonditions. 

E Tide 1S of the mixed type with a period of 24,83 hours 
and maximum range of about 1.89 meters. Field measurements 
саге that tidal variations become larger in the bay 
increasing to 1.98 meters at the center portion of the bay 
ME. 1 meters near the southern end. Currents have a 
maximum speed of about 0.7 meters/second in Zuiga Channel. 

Ae anne as maintained for navigation starting at 
Esa Channel and continuing throughout San Diego Bay. 

The channel depth is 12.8 meters below MLLW from the ocean 
πε Vicinity of the South Bay gage. Other channels of 


Minor importance exist in the southern portion of the bay. 
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PURPOSE OF THE U.S. CORPS OP ENGINEERS MODEL STUDY 

ος Που οἱ shape of san Diego Bay and its only entrance 
Acne northern end, its use as a Naval Base, and the in- 
fluence of a large community surrounding the Bay, have 
EDL into consideration the construction of a second 
Navigation entrance somewhere along Silver Strand Beach. 
NIS second entrance is expected to “case aetredtstime whip 
m ос throughout Zuniga Channel and increasing flushing 
pate of the Bay, possibly improving the environmental 
Memeiv ions Within the bay. A major consideration in eval- 
uating the overall benefits of a second entrance depend 
Meavily on assessing the extent to which it would increase 
mec Tlushing rate. These effects could not be computed 
reliably by the available analytical methods [10], soa 
ШО сс hydraulic model of the bay, in which several 
n Led locations for a second entrance could be tested, 
mes constructed to investigate the effects of each proposed 


Mecation, in detail. 


tee THE HYDRAULIC MODEL ΐ 
The dominant forces in rivers, narrow estuaries and 
Some lakes, are pressure gradient forces driving circula- 
tions which are opposed by friction [12]. In large scales 
of circulations, the opposing forces are combinations of 
10η and Coriolis force, or Coriolis force alone. 
Similarity can be considered as being graded from 


Beometrical through Kinematic to dynamic in the attempt to 
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DC e existing conditions in a scale model. A 
model has geometric similarity if the ratios of all homol- 
Erous dimensions are equal. The model is said to have kine- 
metie similarity 11 the paths and patterns of motion are 
Ewmetrjeally similar to those of homologous occurrences 
he prototype, ana the ratios of homologous velocities 
Ares equel at all times. Dynamic similarity is achieved 
wben the ratios of homologous masses and forces affecting 
Mevion are equal at all times. Complete dynamical similar- 
Mees difficult to achieve in fluid models and, kinematic 
Eas ur:ty Ol necessity, involves elements of dynamic 

en larity. 


ου σα equation of motion can be written in the 


Lorm: 
~ Ξ -p20xV - Vp + okg se 
here: 

ΙΕ Torce 

o24xV Coriolis Ío orce 

ур ο ο ους гаалтепы Гогсе 
oKg μα πο force 

F = fr el io force 


Each term of the equations of motion represents the 


forces present in a prototype and what is to be represented 


1 





mama mocel. Since the totality of forces cannot be easily 
E odueed in hydrodynamic models, a few important terms 
are usually singled out. 

In sealing a model, each quantity and dimension of those 
er cular Terms of the equation of motion used to describe 
pee LOW must be ομαδα Mames model scale. Dimensionless 
I IO5 of terms are used to achieve correct model scales. 
Me numerical value of the dimensional ratio must be the 
Same for the model as for the prototype for simiarity to 
erevall. 

Several dimensionless ratios of terms are common to 
oceanographic modeling, and each one is used depending on 
Mee naracteristics of the flow in the prototype. The 
ice Humber 15 generally used in models of estuaries 
Eu se it involve the ratio of inertial to gravity forces. 
EN "ocprces sre. With friction, the dominant forces. The 
e Lion of flow is primarily down the gradient of 
Pressure. 


The Froude number is defined as: 


Pr = mercial term u B 
р = -—_ چ‎ 
puc ional term 1/2 
(Ly/p) 
Ши ге: 
xs NENAS tac velocity 
U = Characteristic length 
y = specific weight 
p = fluid density 
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Sa pre lilmimnary step 1n sealing a phenomenon down to 
EN ! size, ic is useful tO construct a ratio of units which 
Convention 15: 
A 
E 


A 
P 


R = 


where Av is a number representing the dimensions of a cer- 
mn property of the model and E. the dimensions of the 
Eworlosous property in the prototype. 

"πο length ratio is fixed first in order to match the 
available space or equipment. The numerical value of each 
Gemensionless group constant can be held altering the ratios 
of the other units included in a dimensionless group. By 
... Procedures, time and other unit dimensions of other 
related properties may have to be altered to permit the 
IMC. motion in the model to be regarded as behaving in 
E Same physical manner as that in the prototype. | 

ШО КОЛО О of inertial (velocity and length) and gravity 
is of concern for the Froude Number. Shallow-water gravity 
Waves, such as tides, propagate at a speed OC 54η. 

The velocity ratio R? is satisfied using the length 


ratio (Ry), where 
C 


lA 
-.m./ mm . far i: [RO 
τ ς Is BUD an 


p 


E9 





which assumes that the gravity ratio Sec) is one. Therefore, 
imme length ratio and velocity ratio are left as adjustable 


poprameters 


It-can-be-eseen from this equation that R. shortens with 


mucpPeasing БВ, and lengthens with increasing R 


h J. 

A distinction between the vertical length (h) and the 
horizontal length (L) is often made. If the Froude model 
N:Nocometrically similar to its prototype, Ry = Һу, 

The vertical scale (н) generally must be chosen so 
pvc the least depth of water in which the flow is to be 
εὐ js approximately one em. to prevent capillary 
EE ccs fron being important. Then, it may be necessary 
Dart irom Similarity by distorting vertical dimensions 
euch that R. > E One is led to a choice of scales by a 
p which touches first the minimum depth requirement and 
then an adjustment of the R. and R. ЕЕЕ Пас вае пет сяле 
us mes constant, leaving the others to be adjusted. 

The R, (Aine ratio) is usually the most flexible 
parameter available and determines the frequencies of the 
mulated tides. 

ИШ Мег СО Б Of wind stress on the circulation in estu- 


aries and bays are often important and sometimes dominate 


the tide. To scale it in a model requires more knowledge 
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Si che dynamics of momentum transfer across the air-sea 
mieeriace than iS now available. Empirical procedures 
που be made to suffice. 

ASUS oO ps ol engineers hydraulic model, the 
ton ine characteristics and scales were used: 


mea: 2090. sq. km: 


Horizontal scale (Rr) = 1:500 
Vertical scale (А, ) = 1:100 
Velocity scale (R,) = quu 

Time scale (R, ) 50 
Discharge scale (Ενα) = 1:500 000 


Volume scale ο 2252000000 


Physical dimensions = 35 x 40 meters 


ΠΕΡΙ OF THE HYDRAULIC MODEL 

Mito pero: i low being studied determines the predomin- 
meet oreces to consider for dynamic similarity. The San 
inso bay model used a Froude number scaling criteria. 
EM Pr determining the dynamic Scaling, a hydraulic model 
memeonstructed tO match the scaled geometry of the prototype. 
EN г reproducing the geometry as closely as possible, kine- 
ac similaprity must be obtained. This is accomplished 
by matching the hydraulic head and flow ns at various 
locations in the model over the tidal cycle. This matching 
is generally accomplished by using roughness elements which 


ENEMC 0:955 0f concrete blocks or metal strips. The rough- 


a S elements have the effect of increasing turbulence and 
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ШОО о mean flow patterns. This is a rather laborious 
u (o procedure and usually is the most time consuming 
meet Of hydraulic model building. After the kinematic 
СЕЕ apr3jty is satisfactorily accomplished, the model is 
Presumed to be calibrated. 

Uee eS CIne Conditions were made under carefully 
E Dolled conditions of tides, currents, and simulated 
poBtion input [i1]. The results of these calibration 
ucc called Base Tests, were used to evaluate the effects 
s avigation openings. Thus, any differences noted during 
the tests of the proposed second entrances were attributed 
EU »ences of the plan being tested and not to errors 
Ehe construction of the model. 

ИО ООЛ ОСОО ОТОП on the day of observations can in- 
troduce large modifications to the normal surface and cur- 
rent characteristics in a bay such as San Diego Bay.  Dif- 
ferences from day to day in the values of a mixed tide make 
Ene-day observation a non-periodic function. Therefore, 
Meertications of the field measurements were introduced in 
Mer to make them periodic and comparable to a non-wind 
EN nation. 

Micewcitects Of the souch entrance on tides throughout 
the bay are shown by comparative tide curves for measured, 
base and planned πο. ea Те зесопа 


vce caused a reduction in the tide range at the southern 
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end of the bay by lowering high vater and raising low water. 
Po the south bay gage, the reduction in tide range was about 
NS cm. 

The effects of the second entrance on mid-depth current 
velocity are shown in figures 2 and 3. The opening reduced 
Ehe Maximum mid—depth current velocities at all stations on 
Ens through "D" by 15 to 75 em./sec. The flushing of 
Me” southern portion of the bay became governed by the new 
EN nuce and there seemed to be little interchange of water 
Mer meen the northern and southern portion of the bay with 
Marea 1n the vicinity of Pier No. 2 gauge having very 
Eu currents. Under those new conditions, the entire bay 
appears to have two completely different circulation 
Ecems driven by tides at Zuniga Channel and at the second 
K ce, The flushing appears to be increased in the 
mu nern portion by including à second entrance. 

пае расгсогіал representations of the flow in the bay 
Mager both conditions [1] are shown in figures 29, 30, 36 
EN Tn these figures, the flow is represented at 
EM rcnt stages of the tide, and its intensity is scaled 


the length of the lines. 
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ІІІ. THE HANSEN HYDRODYNAMICAL MODEL 


rhe assumptions of the numerical model are: 


ШО ОШ О e ioo r ous and incompressible. 


Ice rohwsPoSvsbre cand thus, the changes in 
pressure are due solely to changes in water 
n ace elevation. 


K πεεζύτοη berms are vpnorpred. 


Te fluid is in hydsestaticweguilibsium in the 
οσοι σοι οποιο, 


... ΠΕ geographical and vertical variations of the 
ο πο τε σος ους πες Ιςσοῦςςα. 


Applying these assumptions to the equations of conserva- 
irom Of momentum and mass, basic equations are obtained for 
the single-layer model developed by Hansen [6]. Horizontal 


memencvum equations are integrated over depth to give 


Ju 2 an 2 b 

pt fv = - E x г Е ЗСО τίν) 
ΟΥ 9 η 2 b 

A fu E 3m V το τι 
an 3 0 = 

SE + E (Hu) + gy “НУ? 0% 


The wind stress components are represented by 
S rA. 
X X y 


2 2 
W W + W 
ΠΩ͂ y 


+(x) 


enl 


|» 


tly) 


e 





Aa Ne οσο ο Stress components are represented by 


ie? + E and 


vu? + y 


τὸ (κ) Ξ 


LA 


т?(у) = 


Ea 


The various terms of the equations are defined as: 


K = coefficient of horizontal kinematic eddy 
viscosity 


\ = drag coefficient 
Wy ) = wind speed component 
n = surface elevation 
H = total depth = h + n 
х,у = Space coordinates 
UV =- velocity components 
e aeoe leration of gravity 


coto ls parameter 


9 9 
zu ο. 
ð X dy 


Both the wind stress and bottom stress equations are 
Expurically developed. The wind stress terms, as used in 
the model, are assumed to be a quadratic expression in wind 
speed, where A is the wind drag coefficient with a typical 
Value of 0.65. The bottom stress is assumed to be non- 
linearly dependent on u and v, and its formulation, like 
the wind stress, have been derived by empirical means. It 
was originally formulated for shallow water applicetion. 
Over deep water, its value becomes small and its 


applicability is questionable. 
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ООЛО О ЛОТ Села difference scheme is used- for 
e time dependent solutions of the equations. At 
SS boundaries, the values of n, u and v are taken at the 
S1 points rather than from the surrounding points as 
shown in Figure 4. 

υ  τηιϊνίηρ Forces are tides input at the open boun- 
es. and wind at the surface over the entire grid. The 
tide values are computed at each time step and introduced 
ev values of n at each point of the boundaries. 

The finite approximations are given below. The water 
ασ elevations are first calculated using the 


foemservation of mass equation. 


AT AG A 
Mo πι PE £ t 
n imc (n,m) - Ed (n,m)U’(n,m) 


- H  (n,m-1)U" (n,m-1) + B." (n-1,m) V" (n-1,m) 


-- HY (n,m)V" (n,m)} 


The horizontal and vertical velocity components are then 
determined from respective horizontal and vertical momentum 


equations. 
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. CAU = ο. = 
u**At (nim) = (1 - [at рини? ^E (nm) MUT Qum) £V." (n,m)JU" (n,m) 
¿AL 
x 
+ at fv °(nym) - S828, ? (n,me1) 
B 
Ë „2172 (п,п)) + At FE am) 
ZU КИРЕ с uuu c Mae > 
ΔΕ mr [55 ES St 
V (n,m) » (1 - [^t r/Hv Cosme ly V CRM) FU Cam je V Krom) 
ciat 
x+ 2 
= At f U V (nm) _ 26/28 4 ° (n,m) 
t^t 
~ n (n+1,m)) + At Y (Man) 
where 
U” (n,m) = aU” (n,m) t =” (U* (n-1,m)4U" (n*1 m) &U* (n,m*1) 
£U* (n,m-1)) , 
V* (n,m) = av” (n,m) + =” tv (n-1,m)+v" (nt1 ,m)+V" (n,m41) 
£V (n,m-1)) , 
z 
U U (nm) = Т CU* (n ,m-1) &U* (n1 ,m-1)4U* (n m) *U  (n41 ,m) ) , and 
x 7 
V U (nm) = y CV* (n ,m-1)4V" (n&1 ,m-1) *V (n m) 4 V* (n41,m) ) 


im factor a Can be interpreted as related to the horizontal 
kinematic viscosity parameter. In the program it is treated 
as a tuning parameter and is related to the eddy viscosity 


coefficient by: 
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he total depth of water is calculated by: 


ΕΕΔΕ i tti + 

Hu (n,m) = hu(n,m) + 5 {n (n,m) * n ο. mt1)) 
M 1 2 ty 

Hv (пош) = hv tae) + 5 (n (n,m) * n CR m+1)) 


Emu cetfects of wind are computed by: 


К sf at at PULO or зр. 
τ „‚ 1 0 and 
ONE CX 
ΠΟΠ + (W 5° 3P 
аа o 
р ду 


where D. ο ο cE E comence pressure. 55 gradient is 
assumed to be zero for small areas and normal conditions. 
mer stability of the model is governed by the Courant- 
femedrichs—-Levy criterion which says that the maximum length 
of the time step is determined by the grid size and maximum 


depth in the area, 


i 
28H 


AC < 
max 


s particular grid, a At (half-time step) of 5 seconds 


was adopted. 


p 


NN NcoMputatclonesdt area contains small sections of 
greater depths [4], a "false bottom" can sometimes be 
assumed in these areas. (E.g., areas deeper than 500 
meters can be assumed to be 500 meters deep.) This will 
aten result in a considerable increase in the time step, 
Sema decrease of the grid length if total computation time 
mena critical factor. Experiments have shown that the 
ο introduced with the above procedure is acceptable 
some cases for practical applications of the model. 

ο ει το ος «οπου. of three different sets of grid 
points shown in Figure 4-a; the water elevation (z), the 
velocity component and the v-velocity component. Each 
E hese three points have the same coordinate designation 
(n,m). The coastline must pass through u and v points and 
mu Through Z points and the values of depths at any of 
Frese particular points are read in the program as HTU 
and HTV. 

The HTZ points act as water-coast-land designators 
Meine values of 1, -1, and O respectively. Тһе НТО апа 
HTV show the depth values (in cm.), the coastline as -1, 
Ex :nd as 0, depending on their location in the grid 
[Figure 4-b]. - 

che input open boundaries, HTZ points have values 
of -2 (and -3) for the first and second open boundaries 
(if applicable). Outside of the input boundaries, values 


5ο for HTZ are prescribed. 
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Ihe"coeclficient of horizontal kinematic eddy viscosity 
mee interpreted as related to the values of the u and v 
Event components. Hansen states that the computations 
are always stable for values of Beta (l-a) larger than 
E ccand the normal value used is 0.01. However, this 
weent can also be used as a tuning factor. If a 
Meer Value is used, the current speed is generally 
@ecreased. 

In areas where the ο puEdMcUPIDUULOM is irregular, 
meeelcravions and surface irregularities appear in the 
model. This abnormality can be solved by means of a proper 
smoothing of the bottom or the sea surface. This is accom- 
Moneda in Chis model by the smoothing of the sea surface 
EM Scion which is a numerical artifact of the model. This 
E: meme со insure numerical convergence. 

The coefficient of horizontal kinematic eddy viscosity 


Ne tinite solution, is represented by the relation 


EO) 
HAt 


he values of AL(100 m), At(5 sec.) of this specific 

model, and alpha (0.992) recommended by Laevastu for estu- 

aries (personal communication), K becomes 0.9 x 10? em*/sec. 
Bowden [2] suggests for estuaries with a tidal current 


EDU cude of V and a depth of h, 
Ка 0.15 Vi 
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This relation gives a value of K = 0.25 x 10? em“/sec. for 
ШШ northern portion of San Diego Bay. This corresponds 
moral aû value of 0.995. For the southern part, this equa- 
oh gives a value of K = 0.3 x 10^ with a Correspondine 
value of a = 0.999. 

fie use of these high values for a must be done with 
L on because the low smoothing can allow undesirable 
@eertlations in the model. If a too low value of a is used 
Boarder to obtain good smoothing, the resolution of current 
meula be affected. 


As was pointed out in the assumptions of the Hansen 
D 


Bee, the advection terms, u. RS bave been neglected. 


EE 
If a value of 50 em/sec. (about one knot) is assumed for 
Doch u, and Usa Toni cne distance between grid points is 
200 m, then, the advection terms are of the order of ome 


cm/sec. A representative value of the local velocity change 


3% 
at? 


EN Dnoups during which time the velocity changes from a 


ο ο (ОООО СШ ср а At ot a quarter tidal cycle 


maximum represented by 50 cm/sec. to zero. This term then 
EN che same order of magnitude as the advection terms. 
The coefficient of horizontal eddy viscosity (K) is deduced 
fom the horizontal advection terms; then, in the model, 
what is partially done by the tuning process is accounted 


mer Tor the values of the neglected advection terms. 


E SELECTION OF THE GRID 
The selection of the grid size is usually based on 


meqduarements of details and accuracy, and availability 
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of computer core memory and time. For open areas and round 
shaped smooth bays, where the expected velocities and direc- 
telas are smooth, it is not necessary to use a fine 
Mattei areas where topography is of primary importance 

EE une mesh is necessary for reliable results. 

The flow in San Diego Bay is governed by the topography 
Me area; therefore, a small grid size must be intro- 
duced. To cover the entire area, a large array would be 
necessary which would require a large amount of core and 
amputation time. 

VM the other hand, due to deficiencies in the boundary 
KS l Ons oi the model, inaccuracies of the computations 
EM uno near land grid points; therefore, if acceptable 
meters are desirable in narrow channels, the grid must 
meevide enough grid points across these areas. 

w al erid or 100 x 100 array was lain over San Diego 
Bay. Because of the difficulty of handling by the computer, 
rea was divided in two regions with sufficient overlap 
Ensure proper calibration. 

The way the boundary conditions are set at the open 
boundaries produce errors which propagate throughout the 
EM cent grid points; therefore, the overlap area must be 
oí considerable size in order that the match section be 
distant enough from both open boundaries that these in- 


ences are negligible. For this specific model, a long 
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Narrow portion of the bay in front of Coronado was 
Hosen, assuming that jl Ti αρθρα flow is parallel to the 
ss (L line and free of big, eddies. 

MO is ton glves two prids: a northern grid of 
53 x 30 designated as grid (1), and a second southern grid 
of 58 x 40 and designated as grid (2) [Figure 1]. 

Me Cond portion of the bay was run under two differ- 
conditions: first under existing conditions (model (2)) 
сие secondly, with a second open boundary at Crown Cove 
(model (3)). The division of the bay into two different 
Mete ls will not allow for computation of velocity at the 
memonern entrance for the case of the model (3), but it 
ο лүе а food indication of what the future conditions 
will be in the northern channel as the point C-2 (see Fig. 1) 


[meee Computed in the second and third models. 


Peel i DE INPUT IN THE MODEL (1) 

In the northern model referred to as model (1) (the one 
with two open boundaries), tidal values at both openings 
Baszrırst introduced, The time inaccuracies of the proto- 
type measurements forced water in and out of the bay at 
m@emewo ends simultaneously. This procedure causes a non- 
MNR Dian state in the bay with corresponding erroneous 
answers. 

A new procedure was tried by leaving the second boun- 
Ev open without a driving force, letting the water flow 


ana out according to the propagation of the tidal wave 
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τ. ΙΕ τΙτ5ς boundary, Continuation of the surface of 
meenvaver and currents for each point in the free open 
mary was prescribed to insure the flow oí water through- 
is end. For the values of height, this continuation 
was obtained by leveling the edges using the concept of 


a “lubricated wall", 


Z(NEH,M) = Z(NEH-1,M) 


Ие рсе to continue the surface slope resulted in distur- 


Rumes being generated in the model, 


Z(NEH,M) * Z(NEH-1,M)* [Z(NEH-1,M) -Z(NEH-2 ,M) ]. 


νι πρ the edges and letting the model compute the values 
Of Currents itself resulted in too severe a decay of the 
[memes Of Velocity at the second boundary; this resulted 

in no flow through the free open boundary. 

ЕОР спе соггесс пон ог water throughout the free 
open boundary, a continuation of current (u and v) was 
Eucserjibed with the assumption that calibration of the model 
ὃς Obtained modifying the values of this slope by 
means of applying a smoothing constant (J) to its value 


Meeeure 5], 


V(NEH,M) = V(NEH-1,M)+J[V(NEH-1 ,M)-V(NEH-2,M) ]. 
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Warliations of the values of the currents in the entire model 
were obtained, but great disturbances in the water surface 
were introduced up to a point that back-flux was obtained. 
Figures 6 and 7 show different results obtained for differ- 
ΙΙ 11185 οἱ α 1η the continuation of current used at the 
free open boundary. 
These disturbances in the surface can be reduced by 

EGER ue an inverse correction to the slope of the surface 


at the edge 


Z(NEH,M) = Z(NEH-1,M)-C(Z(NEH-1,M)-Z(NEH-2,M)). 


ο er calibration was attempted because this procedure 
Banner realistic, but a proper calibration may give a good 
mt with savings in computation time. 

Because no satisfactory answer was obtained from a free 
open boundary, the idea was abandoned and the addition of 


momappended pseudo-bay was tried. 


C. ADDITION OF AN APPENDED PSEUDO-BAY 

Me addition ol an appended area at the end of the 
Ml (1) creates- an additional tidal prism which will 
Bee the volume of water passing throughout the first 
model to increase. This new area, Figures 8 and 9, was 
obtained by filling the empty spaces of the grid and con- 
Maine them to the end of the model. This procedure does 
NS DcPease the core in the program, but increases the 


Pine Of Computation by about 30% for this specific case. 
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Areas corresponding to 82% ana 100% of the remaining 
area of the bay were tried. The results, nau coc. ls 
aná 12 show an increase of water transport and current values 
with no appreciable modifications to the values of the 
water surface. The increase in volume το σος ατα СИТЕ 
rents were proportional to the increased areas. The values 
of height remain stable for uc ditlorent cases 

Different values of R (0.003, 0.0028 and 0.002) were 
applied to the model. The results show that the model is 
not very sensitive to the variations of this parameter. 
After comparing the results obtained at the points A-2, 
mee and C-2 for each of the values of R, 0.003 was adopted 
because it seemed to better fit the measurements at these 
points in the prototype. The discrepancies at A-2 were 
attributed to the inaccuracies of the Hansen numerical 
model near the open boundaries. The results of the model 
(1) in current, water height and net volume transport are 
Myewn in figs. 13, 14, 15. 

Bepictorial description of the flow throughout the 
entire model (1), with the appended areas, is shown in 
eure 9. The real part of this model and the times of 
0., 3., and 9 hours are shown in Figures if Por 
The direction of the flow is indicated by the arrows, and 


its magnitude is scaled by the length. 
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D. MODEL (2) 

From the results of model (1), tidal values were ob- 
tained for a point at the middle of theschannel in the 
vicinity of Navy Pier Tidal Station. Theses elues vere 
Produced as input to model (2). The computations were 
attempted with the same values of R (0.003) and a (0.992) 
as in model (1). 

Using these prescribed conditions, oscillations in the 
current and tidal values developed (see Figures Io oand COE 
These oscillations appear to be in phase and their magni- 
tudes increase from zero at the open boundary to a maximum 
value of about 30 em/sec. and 15 em/sec. respectively at 
the end of the deeper channel at the near end of the bay; 
and from there, the values of the current reduce proceeding 
to the shallower and southern part of the bay. 

The oscillations seem to be activated as the tide 
approaches the LLW, increase their values throughout the 
next HW and die out by the next LW, repeating the eyele 12 
hours later when LLW is reached again. The period of the 
oscillations is greater than 2 hours. 

The causes of the oscillations could be some kind of 
seiching manifestations in the model, ue water tiae 
or oscillations characteristic of the model produced by 
attempting to work in extremely shallow water without 
introducing special prescribed conditions for these small 


К лез [8]. 
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impossible mode of seiching can be calculated assuming 
that model (2) is an independent bay with a length of 14. km 
Eus medium depth of 4. meters. This bay will oscillate 
mL fundamental period of about 2.3 hr. which agrees with 
m resulting values. A Fourier analysis was made of the 
tidal input of the model (2) and the calculated water height 
sat South Bay Gage. Calculation of the Energy Density Spec- 
trum was performed to compare both and see if energy from 
some of the harmonics from the tidal input have excited 
the oscillations in the bay. The spectrum shows (see Fig. 
23-1) that the oscillations have periods between 2.06 and 
Seo (with a maximum at 2.3) hours and that the tidal input 
negligible energy in these periods. This seiching may 
be activated at the LLW and damped out by the smoothing 
parameter a. ‘The existence of such anomalies are not 
shown in the published records of the observations and its 
EN nee depends, in case of its existence, on the time 
F εἶπςσς of the data obtained. 

ШОМ Оо рл Ту ОГ а shallow water tide seems to be 
Eu conable because in each tidal cycle, the oscillations 
TN TO be activated by the tide when the depth of the 
rer is at its minimum value. This shallow water tide 
produces a high frequency harmonic [3] which can be ampli- 
fied by the model and damped out later by a. 

Because the model has linear terms only (except in the 
meevom and wind friction terms), it is not possible to 


miler energy from one harmonic to another; then, the 
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amplification of a high harmonie from the shallow water 

pde is questionable. Ое а5 0 possibility ean bee defi 
πεν of the difference scheme. As the tide is a continuous 
NET ype propagation, errors are introduced in the time 
differential computation that could produce a harmonic that 
could be amplified by the same model. 

An.a of 0.99 was introduced in an attempt to eliminate 
um oscillations by increasing the smoothing, but no 
improvement was obtained. Lower values were not introduced 
Because larger smoothing may affect other parameters in 
Ene model. 

meecemeca ior a iow value of Alpha in the model to damp 
Ku L ese oscillations makes the previous estimation of 
ШОО of a not applicable for the southern portion of San 
Diego Bay. In the model 2 and 3, a minimum depth of 4. ft 
memesavablished to prevent that portion of the bay from 
πρ ανν at LLW. This minimum depth was increased up 
EE it. to test the probability of shallow water tide 
[eon the result that similar oscillations appeared but with 
larger amplitudes, (see fig. 21). This new test made the 
possibility of seiching the most feasible mechanism causing 
mae Oscillations. 

Model (3) was finally run with an a of 0.992 and an 
mone) 0.003 and the results of current and height of water 
οπού in Figures 22 and 23. A pictorial description 


Еспе Circulation of the Model (2) at the hours 3 and 9 


D 





in cen/sec 


Current 


Height ín cn 











А; Ñ. : | | 2 ра lunar 
oe. PES X Р T $ > 
x ag ° 
ok ος 
— s 
` 4 Numerical 
| Minimum depth 
E^ 4. Ít. 
12, ft, ----- 
CURRENT AT С-2 
MODEL (2) 
12), 
90, 
60, 
30, 
m 24 
Lunar 
-30, τ time 
es 60 ° K / 
e 
- 90 е iau 


HEIGHT AT PIFR N, 2 
MODEL (2) 


FIGURE 21 


60 








Prototuse ------ 








«90, Hyde s EL d 
Numerical A 4 
CURRENT AT C-2 
MODEL (2) 
120, 
90. 
о 60, 
0 
Twn 
ç: 
Ὁ 
E 
ec 
Ë a 
L 
т CAE 
fa 
з 
O 
«60, MA -- 
Hydraulic 
-90. base ------ ----- 
| CURRENT AT D-2 Numerical 4 ——4 
MODEL (2) 
FIGURE 22 


61 





in cn 


Height 


in cn 


Height 


120. 


0 


60, 


30. 





“60. Proton ------ 
Hydraulic 


Ваве - 


Numerical 4—4 





«90, 


HEICHT AT PIER N, 2 
MODEL (2) 


120. 


O. - 


timo 


N А Hydraulic 
Wa / Date „a 
Numerical 4—â 
HEIGHT AT SOUTH BAY 
MODEL (2) 


FIGURE 23 


62 





E Shown in figures 24 and 25, and of the entire bay in 
which model (1) and model (2) have been matched is shown 
in Figures 26, 27, 28. In all of these figures, the 

=m Lion and magnitude of the flow are represented by the 


direction and scaled length of? the arrows. 


E. MODEL (3) 

Several alternatives of a second errante anda Mfes Non 
Meme tested by means of the hydraulic model conducted by 
mem, >. Corps of Engineers. Because the purpose of this 
work is mainly to make a comparison between different models, 
only one of the alternate second entrances was examined in 
Bacenumerical model, and no diffusion study was made. 

Model (3), with the proposed second open boundary at 
Crown Cove, was run in a similar manner to model (2). The 
maemo Lor the northern open boundary was obtained from the 
results of the model (1). The same oceanic tidal values 
cured at Ballast Point, in Zuniga channel, were chosen 
ne input at the proposed second channel entrance. 

The same currents and tidal values developed in the 
model (2) using the preseribed conditions with a and R of 
0.992 ana 0.003 respectively were used. Their character- 
τος and interpretations were E E the previous 
Eon. No other values of the tuning parameters were 
ved because the-purpose of this last model is to test 


model (2) with an additional open entrance. 


63 





< ғ э > ? 


< Σ > FP = æ — = 3 


Ф Э эу „> э» - m » o ә 


* € r ο ος ο» 
< < ДИЛ ee m 
το. τς ο E 

RI E A T, 


€ < 2 > эъ — = 


< - A IA A A TEE RE ETS 


< DEA O O کک ص‎ ы U U UU I; >> IIA э ж > сс 


є Оу A E um ; =m ene eds ЕВ O m= m د‎ eee » z = 


€ WW M ia N RE p ج‎ < < = 

€ UE ERE — —  — — _ > د > > و‎ > 2 2 < = 
A η ا ج سے ج ج — صد سے سے‎ EE e œ 
TE νο ο οι. ae см 


с e 
E ГЕ Ur er E E UR . e 
А К τον ЕВИ Е ТАВИН ου. t © u ef - 
Z 8 е 


| 
| 
i 
| 
| 
| 
v 
N 
ч 


| 
| 
| 
| 
| 
\ 


- Ὥ — — — - 





A A A A A 


t A A A A A 


CURRENT DISTRIBUTION IN MODEL (2) 


Зеб 


TINE 


FIGURE 24 


64 





€ < < c CE < 


¿< “ - 

ae 

е сг 
CUI c DE 
EG CE TT 
AE ASS 
` > IO A ως 


= єс e sc ε <i < 5 





we ee m n a m A ~ =. 
О Е c AAA ARS EA AAA 


TA AA AA 
а 


< 





E a к — 77 «t ac 7 "— — 
S A A A ος gee 
« A: AE —— u T Ta mm سے سے‎ 
< προ σα e с = = = 9 а œ 
τ ( ж A Уа = € - - 
< — ς᾽, 'α-- — — «7 a7 - 
< а а — — e 
< 
‹ o a АЕ СР СС 


< 
nr 
e€ N 
N 
T 
/ 
a 
А 
/ 


c¢ > œ >o a œ œ 


= 9hr. 


CURRENT DISTRIBUTION IN MODEL (2) 
AINE 


FIGURE 25 


65 





Ф «а . - я e s 5 * * * 


eevee 
° a - = . 
wi mu m f£ € 
= e pr --- - - r p e 
ы =f lar © = = a a a s 
я .- rn. « A “° $ 
<“ e + 
4 
W. ο ο 8 4 . © = 0 - ` “ = 2 А Κ΄ a a 6 e e 
£ τα La 
* = a r a P --- -—— mee we = = a 
>= 2 s чь. m = Ф 
PE >” ~. ar ao — E ےت‎ == = == -- 
- 
~ 4 
- pu “ч - ——-——— — — -- —— م‎ ΗΟ, ч > w un 
а -- e y e 
- a 

= ت سے سے سے - - > سے = pe‏ - کہ سے عور —- - - - -— 

акы a a CA ο αμ m er “Y E . 
ET — „чыла σας E LS í T: nn 
Kun -- —— w > Os ہے‎ ee νε 
. 


a 
a 
« 
кез се s ear aea а Жо к" 
£ 


< 
ee 
= > € -— [u € M — p s= = я” ر‎ ΄ 
Z < a <“ “ < - 
. w AA A A ar ca ar am at Um ^ 
a pens - 4 
«К  « а 4 : - æ 
те - — a —— oe eg ar cr up μι -- РА 
e P a Te a - 
------ᾱ---------------α"'.ε΄”" α΄ α΄ °a = 
e РЕ 
SSS m cM A A οι ο 
y un = νε “> ο =“ =” = 
απ να τα . ο 
e Pane Oe ee α΄” = = ہے‎ 
2 А ΄ 


κ΄ 


ο - 9‏ ات ےج مدصت د کے ر ےو 


* Ihn ο ο Der De “weg - 
| / y O “< ν. ῳ че 
ft fe ў οκ πμ. τας. == 

= ατα p ee e = = 
Effy i ç ` ` 
* А / 7 al Е r = nm u. 
FEAT A EF ا‎ 
IDA [er x ΝΝ-- Sn 
M 4 4 / j 
E eee, pi al ë © © t ` eee 
Lok ye yer a 
. o. MAA ^ f т + - . 
/ 74 7 4 са с € i 
- > £ í o ο es eee e а 
رر‎ ae “ “< 
е Z Z lA v rd a 
ay " 
- ae . . 9 9 « ο е 
zZ 2 
- ^ у ` . . . 
“ 
. ΄ . 9 
eA. y 
wer . * 4 Li . 
/ 


N 
- 
` 

я я * * 


"ut A EA pf if) , 

т d سی ار سے لے م‎ F е, . 
TN P 

} Jo» η z 
4 ` 


ج -  —_‏ ——— — سم س ج سے سے 


CURRENT DISTRIBUTION IN MODELS (1) and (2) 


= 9 hr. 


TIME 


FIGURE 26 


66 





“m Q > мә 


< < — 9— 9— مھ ص 


A A а а i e 


po -.- 


- ч» - y 


AO 
— 
- 


» 


2, 
7 


-- 


— — — — a 


e‏ د 





ee 


و 


f 


т 
⁄ e Er tmt 
AAA 
⁄ Z “z —.“.2 2 
EAT е 

^ ^ 


JA 


° 
. ç 

< o 2 ο 

"T. 

- ` - 9 a 
- = LA ЛУ. 
=> = = 2? >> ” P? 
= = = 7 P ” 2 
en» = # » > 


= = =“ “ > ” 2 
- = o» .. 3 
- 

очо оч 


το o, © е = LJ - 
. - «a 
° e = = Ó a e - » 
s e ewe =“ <“ °“ “< <= 
. = = = © =Ыы o m ~ m a mn < T“ = 
. T ж ж = ® = ы = ж m“ > ж = m” m” <= 
э гюл ο ο a a lm s E s m оф = 
` 7 ο ο ν E EEE CO == ” #7 
`. ο ο μμ с ert esa as es em # - 
. DET AM MED D Nc p ee жое = ш = = = 
+ 
- A Жш т ss е = = = ж æ n 
# 2 7 
e 2 WP mL ο ο. eo. ”»o = 
ΑΡ 
. ? NSEC ECS 
A 
- ον. | — а чш ш э Tü =“ m = э © . o 
Jy e 
* 2 SAS Ie . 
a FUA 99. mn < eo чт т - - 
2 № * 2 
: αμ... LEM - 
E T 
- Γ.Ν. 2 - 
FA WR η S XL proe 
- WE DE d 4 L] 
x . o 
RS» yo x 
` PY ΄ 6 
Gi v е œ ee «ч 
f 
= » πλ N) >. ο e v 5. 
΄ ر‎ 
ا کے‎ a a Re т ο 
2 
a - 


- vw. 
е = w m o = w 
@ w we Few e 

Ñ e ç ? 2. „ 
wen." ee 
= =» = = = == 


RENT DISTRIBUTION IN MODELS (1) and (2) 


C 


ААПТ 


ТІМЕ 


FIGURE 27 


67 





f d a 


r 


Рол л < 


r 


f é ur 


, 


* “< < =< < <q & * * - ο . 
т 4 ш αι a 
` - ow o ow ` г MAA + 
ο а WE nM ME ET S - 
ο ο э» лае» „= ay M 7 « 
eS 13S E ME ي‎ 5: 
LEES IE SE Ww xq TN * 
a. 4 Se а ш оа AAA 1 Ы 
> . .. . ^ е 5€ MAA - 
т а «= - ` N - 
. € - ` τα .. 
. . - Ө ро Ç AU $ ° 
am rN Ne we E 
* < + E E on 
@ < “< < ча 1 
ж 19.76 ` « 
ο 2 ба 
. 
ται « 4 
a WT S TN 
Е e er Re: . «ο 
4 > > v < w < э evt 
К с ОЧЕ . 4 
a р > w w wt 
. £ = $4 44? # =“ m ú w v 
. 5 x S SIE Bb сч « ο : 
. „ = N к, 2 RT 
4 to. ` 
e o a ον ο ¢ QQ Moa ' ` . 
| оо : . . 
t . а у : * 
. . 
2 o Й ο. 
. - - 
«. Y # ¢ : 
“ Z^ . . : 
a >» Ë 
. ¿ f t 
я $% 
. ο í ç £ 
. Mu 
ЮГА . 
« G , 
a s 1 ° 
. . > 
.. ~ ` 
. 
° ¢ 2 
- ^ 3 E 
a > ams pu 
a . ۾‎ Е 
а ^ ν κ a 
A s T . 
е. аа ۴ $ 
9 4 ғ $ 
“ # ' 
.. 2 م‎ 
a ⁄ +£ f 
. , 
* 5 P 
a ΄ 
. r 
. 
а = а ә а e.a 4. Ὁ а 
> 


CURRENT DISTRIBUTION IN MODELS (1) and (2) 
TIME = 15 hr, 


FIGURE 28 


68 





The results obtained by the U. S. Corps of Engineers 
Enn Figures 29 and 30 are compared to the results in 
ENT enbrs and height of water for fis model shown in Figures 
E and 33. Pictorial TES nta tenso the current 
distributions are shown in Figures 34 and 35. The same 
Iu vesentation obtainea from the U. S. Corps en Bogineers 
паис тац1їіс model for the complete Bay with the second open 
boundary are shown in Figures 36 and 37. In both groups, 
mie flow el waver throughout the proposed opertenerance 
ELM Lo be" nesdomineting factor of the circulatzen in the 
Exnern portion of the bay. The northern portion is 
EI used by the northern circulation and a very ii inite 


mest area is located between the two systems. 
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IV. MERITS OF BOTH THE HYDRAULIC AND NUMBERICAL MODELS 


A comparison between both models cannot be made without 
the results of a proper test of the models under the same 
conditions and on an area where ο. measurements of 
actual conditions have been made to compare with their 
results. Тһе merits of each have previously been-enumer- 
rear ana the goal of this work is to evaluate them. 

emo re ordering the merits of the numerical and hydrau- 
models, it is good to recall that for this specific 
case, the numerical model in discussion (Hansen Model), by 
definition, is a barotropic (p = constant) single-layer 
ШОО апа the hydraulic model used a constant density 
fluid. The hydraulic model was based on the Froude simi- 
larity concept (inertial and gravitational forces predom- 
inate). Although the discussion is restricted by these 
Eds tions and are not applicable to general numerical and 
hydraulic models, these are the most commonly applied types 
of models. 

The spatial resolution of the hydraulie model depends 
on the length scale and the measuring capability. In the 
numerical model, topographic features are represented in 
Ms ot grid points. To obtain a good representation of 
May metry, small grid size is needed that results in an 
extremely large array. A compromise between spatial 


Resolution and core size must be reached. 
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mMetrcame resolution in the hydraulic model depends on 
E usunilitude criteria and can be scaled to obtain a re- 
liable result. In the numerical model, it is governed by 
ШО v llabillity of computer time, i.e., the smaller the 
Ben step, the more computer time required. 

ime scaling in the hydraulic model is given by the 
mia ctudescriteria. <Asevas pointed out in the section 
dealing with the hydraulic model, these criteria compromise 
fees Causing distortions in the model and oblige the use 
EM ome Specific dimensionless number for each specific 
case. The introduction of distortions in the model limits 
ЩЕ N5es because some other terms ignored in the selection 
of the dimensionless number may become important. 

in the numerical models, the flow is described by the 
SOT ion of the equations of MOULON Whtem includes terms col 
eME rent orders and degrees. Each term describes some 
11ο characteristic of the flow. The characteristics 
ure flow will depend on the shape, location, external 
forces Ca. Бану тесу rand properties tonales] ud. 
MUN CSseribing these flows by tbe equations of motion, some 
sehe terms become dominant and some others can be consid- 
' 'ΡΤΤΡἼ5ΌΙ1Εε or included in some coefficients. Іп the 
case of the Hansen Model, this was accomplished by neglecting 
E rEnon-linear momentum flux terms. The choosing of the 
proper terms becomes invalid if improper considerations are 
Some. Then, an improper scaling can invalidate these 


selected portions of the equation. 
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Eu" cmores. the finite solution of this equation does 
mere permit the continuous computation of the flow. The 
Ewputation is done by time and spatial steps that may 
Μ᾽ large errors in the model if proper precautions are 
МО en. For example, the coefficient of Eddy Diffusion 
ШИ леа їп the finite solution in terms of spatial and 
me SLeps. Then, its value becomes affected by the grid 
КЕЕ Капа time step. 

ШИЮ ing data irom a hydraulic model is difficult 
ana limited by the instruments and the distortion in the 
scale of the model. Sizes and numbers of the measurement 
instruments are limited, and the region where these instru- 
ments ean perform a measurement in the model is restricted 
nel? sizes. In the numerical models, the character- 
Eu oí the flow is computed at each grid point; thus, 

K available at each point and at each time step. 
КОГО О леб in the differential solutions introduce errors 
EN Lese erid points in the vicinities of the coast lines 
EN uut boundaries. In the hydraulic model, equivalent 
Os exist in the input boundaries; good resolutions have 
been claimed for near shore regions that can be question- 
able because difficulties in the measurements of quantities, 
mW as current speed, in these regions makes it a not 
always tested assumption. Besides, surface tension becomes 
important in these regions that compromises the model. 

The cost and time are some of the larger limitations 


in programs and research. Hydraulic models are expensive 
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Mare months and sometimes years to achieve a satisfactory 
libration. The use of the numerical models generally 
ταδε considerably these limiting factors. A limitation 
Seumerical models is aen Τεστ αν ο o | το 
Hare not always available. 

s (cC älrfusien or dispersion of pollutants has 
meen accomplished in both types of models with comparable 
results [1 and 5]. In the hydraulic models, this is achieved 
ааз гесс measurements of dye materials introduced into 
ас model. The conservation of mass equation is used for 
meso computation in the numerical model. All the previous 
measurement and computation problems are taken into * 


E sderation for a proper evaluation of these factors. 
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V. CONCLUSIONS 


Eacseompar3jsons between the hydraulic and numerical 
models dealing with the modeling of San Diego Bay and its 
paoposed modifications, are considered separately for 8 
Cu Cases and for each of the different tests involved in 
this study. The Hansen Hydro-numerical Model was applied 
une Bay to compare its solution against the U. S. Army 
TS oi Engineers hydraulic model. The description of both 
models and the manner of how these numerical and hydraulic 
models were applied to San Diego Bay nas been discussed in 
detail., 

περνά]. the numerical model is more easily handied 
memmeonpopraphic changes than the hydraulic. The latter 
EE 15e stored ior future uses after being run as is done 
Du whe numerical model because of space limitation. A 
Comparison of time and cost is not possible because no data 
Ia vallable for the model of the U. S. Corps of Engineers. 
ШС Nree numerical models consumed for final runs a total 
2g hours 15 minutes CPU time in the IBM 360/67. Adding 
average of three times this amount for calibration, the 
cost becomes approximately $21,000 (U.S. doll po 
would take about те weeks for a trained programmer to 
Sewn the complete answer for this specific case. 

The added easier possibility of including the input of 


Wind gives to the numerical model better ability to represent 
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the prevailing conditions. The same numerical model can 
änelude the wind. Usually completely different hydraulic 
Meeels must be built to simulate wind because of similitude 
memes дегасіоп5. 

me hydraulic model has better space resolution than the 
mumerical medel. A method was tested for transferring 
arty values from one numerical model to another to give 
r model the ability to divide long estuaries into smaller 
КИШ ЛЫ improving in this way the resolution of the numerical 
model. With the use of this method, two different but de- 
pendent models were run to represent the complete bay. A 
Euctrrebvpilon in the configuration of the bay was chosen for 
EN vsrlapping region. A Similar procedure is used in the 
E xouriec model to save space in long hydraulic models. To 
Eve 3t, the bay is distorted in a narrow location in 
acfold-like procedure. The model is calibrated at both 
ASS of the curve. This procedure gives the possibility 
ОИ ООП with a model of a smaller scale that gives more 
Mets of the area, but special care must be taken so as 
WTO introduce errors in the time scale. A similar pro- 
K сап be done in the numerical model if the excess 
area is small and can be accomodated in the free space of 
νυ grid. The possibility of running two completely 
different hydraulic models seems to be impracticable 
E Se of economies. 

The numerical model reproduced very well the sea level 


cur rrents for the cases tested, provided boundary 
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treatments are correct. The boundary treatments sometimes 
ENS not easily achieved and are a matter of decision in 
meer TO represent properly the coast line and channels in 
the grid. Representation of sea level seems to be equally 
merece in both models, but the current results seem to be 
better in the numerical model except near boundaries and 
maonelines. The seich-like oscillations developed in the 
memerical model must be taken into consideration if the 
Medel 1s to be representative of the total bay. If only 

ο ρου of the bay is under consideration (for example model 
(2)), the sieche oscillations should be smoothed out because 
Mey are not representative of the bay. 

No comparison can be made of net volume transport be- 
cause Me valuation in the hydraulic mođel is available. 
These computations were made in the numerical model based 
on the computed velocities. The computed values seem to be 
E able enough. 

It was pointed out in the section dealing with model (2), 
hat the boundary conditions іп the numerical model must 
be given special treatment when the area becomes dry during 
low tide. A minimum depth of 120 cm was fixed in the numer- 
ml model in order to prevent any problem, and the southern 
ΜΙ. of the bay was eliminated because its depth was 
below 50 cm. No special problem seems to appear in hydraulic 
EX. for the case of portions of the area becoming dry, 
but special care must be taken if areas become shallower 


than 1 cm. because then surface tension is important. This 
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minimum depth had not been discussed in previous litera- 
Mire for numerical models, and no attempt to evaluate these 
Pemoitions was made in this study. 

In the model with the proposed second entrance, the 
114! assumption of increasing the flushing of the bay 
Eu sco be reasonable. The circulation of the southern 
mortion of San Diego Bay will be almost entirely governed 
Ne second entrance. But, from the computed net volume 
ИЕ τ in the section AA2, it is shown that the inter- 
КО ре Of water between the ones southern portion 
DM Hue bay at this section will be almost negligible. 
malig Lhe pictorial representations of both models for 
w pecial case, an area of decrease of current speed can 
КО iced between Pier No. 2 and D-2 gages where both 
models indicate that the flow of the bay will be divided 
into two separate systems and the projected second entrance 
Du provide the water for the southern system. This would 
result in decreased dispersion of pollutants introduced to 


E bay in this region. 
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READING THE CONSTANTS AND DEPTHS 
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RETURN 


400 CALL PLOT(O.,NN/GPPI+3.,-3) 
END 


om, 

= 

=. 

er 

CO 

UZ 

Qw fl 

m>- OM 

> *3x 

ο 5 Cr E ہے اسم‎ 

Kom CN m № >< 

w2 сү e 

uu гү t + | 

a= oT 

XX ~ li eC) 
4 WU) 

ΠῚ O ~ 9 4 3C 3t 

πο ια — --- 

,O>onn > > 

H- $6 t ee — — 


207 0020) 4 Онан |} li = 

OE ΠΠ AD c 
DY A Ly O AO A A Le” 
(O cx O LO t) O >< >= >< >- U 


© 
© 
єч 


(X| Y CO SI) 


) 
) 


XXX TET TD 
ХХ, ҮҮ-ҮҮҮ, 2) 


NN 
HL | >< е е 
а. x + >> 
© XX na 
бү: w >< >< >< 
<Í has mr wa 

м 

ul д LALA 
ZH TE TIC EE 
t> oX% © ϱ Cats 
POSEAN CO SE = 
O LO O OA оа = 
Ot >> í | X 
(ү AAA A TA 
ee ЭС ИПИЕ )Ε-:Ε 


DWI >x >x eK EH Z 
NADA XX>X>ODODOLu 


15] 





"SAN DIEG!» 
š 


"> 


RS 
t,t 


(ΟΝ. Του 
ЕЕЕ 


2 


T 
0 


e 

ei 
zn 
30 
< = 
سے کے‎ 
>+— = 
لے‎ =)٦ ساس‎ 


e» 
e 
en 


WATER НЕ, ІСН 


= je ب‎ e» 


Lt 


d 
A ل‎ 
ANY < 
eos X 2 
ΙΙΙ СУ e 


(2458540,58,CL L1: TITLE 8:8, L TC] 


Ше I) > elf "ANO 


LL? je = 
Link 
mO — 
pew Ht I fF <= 
2U «c 
OZ < 
A ο 
аа 09 «т 
кн ошо 
NOJ e = 
oa 


Exc cru. d e 
NND Yer 
e 9» t iy lli Xu) 
lH NH | =D -OD 
مہ ہے‎ и Ом осоо << 
لح اہ سس س‎ YE HOO ال‎ 22 
οι ο 7 Ang 
Eee IO JOO! SOW 
O IO E O CIS 


O O 
O O 
cN سم‎ 


2,2 





BIBLIOGRAPHY 


U.S. Army Engineer District, Los Angeles (Preliminary) 
Summary Report. San Diego Bay Model Study. Hydraulic 
Model Study. U. S. Army о аши ne ineer waterway sue реті шепте 


Karton. Corps of Enpineers ,  Vielshure In 7 2 impp 
June 1971. 


Boden. nF. “The Mixing Processes In 3 Tidal Estuary," 
J. Air and Water Pollution, Pergamon Press 1963, Vol. 


anam ар art ο Gye cm e ne we re me 


7, рр. 313-356. 


Deoason, A. T., and Warburg, H. D... Admiralty Manual 
ο) Оез, Her Majesty's Stationery Office, 1941. 


Mie Laevastu, T., and habe, K.. A Deserıprıon of the 
EPRE Hydrodynamical-Numerical Model, U. S. Navy 

E ronmental Frediction Research Hhacilaty., Technical 
Ber No, 3-72, March 1972. 


Dr, Taevastu, T., sand CDR Hamileon, G.D,, Compuvation 
er Flushing and ρε i ο Pollutants in the Pearl 

Harbor and in San Diego Bay with Hydrodynamical- 

ИШ rical Models, Us Se Navy Environmental Prediction 


| covon paclisty.- Montepbey. Ca buf ME 

















E Laevastu, T. and stevens, Pas Application or 
Numeprical-Hydrodynamical Model in Ocean Analysis/Forecast. 
ir. The oringle-Layer Models of и Hansen 1556 
Numerical Weather Central, Monterey, Calif., Technical 
Note No. 51, July 1969. 


DAGA Laevastu, Tes and Си аа ао ОЕЕО о Ола 
ea: Current in the Strait of Florida with a Iyarodynam- 
ical-Numerical (HM) Model, U. S. Navy Environmental 
Prediction Research Fac: Facility, μουν, April 
more. 


Bine, Hs G., Investigation о о оле осевееваи 

u ov Water Areas and Estuaries,; Γη; τ᾽ bur 

Meereskunde Der Universitát, Hamburg, Proceedings 

тр. on Coastal Geodesy, Múnich, July 1970, pp. 439- 
De. 


Remming, H. G., Hydrodynamical—Numerical Investigations 


m oeo rizontal Dispersion 0O Serra area Der Ei 
Ú o [IL Fúr Meereskunde Der Universitat Hamburg, Report 
DOS. Feb. 1971, 712 pes. οἱ figs: 





nS 





NO. 


I. 


E, 


ποσο ος οσα Πο μαμα b no ο ο ου ο рсе 
NCU Second Entrance on TeVe liste Scherzereristles 
um її Diego Bay οσο ον PAO Veena Com crence 

τ pjne Pollucdom ono ics Boe cO К лев тосев 
and Fishing, Rome, Italy, 9-18 December 1970. 





E pdpup. d.—Qu.. Johnson, M. Wicsandcbdomings bi. 
The Oceans, Prentice-Hall, Inc. 1942. 


LEM Doc. S.An Inbroduetion tooPhyselicaleoceenoPraphys 


Addison Wesley Publishing Co. inc. 1962. 


ILS 





|. 


л 


ШО. 


ENS ER BU ООС ЕЭ 


Mo h eos 


Defense Documentation Center 2 
Cameron Station 
Mexanoria, Virginia. 22314 


brary, Code, 0212 ο 
Naval Posteraduate School 
Monterey, California 93940 


Department of Oceanography 3 
Naval Postgraduate School 
Monterey, California 93940 


CDR Leopoldo Salas R. _ ; 2 
Direceion de Hidrograiia y Navegacion 
Apartado 6745 Caracas, Venezuela 


Escuela Naval de Venezuela Js 
Escuela de Posgrado 

Meseta de Memo, Catia la Mar 

Venezuela 


Mis Edward E. Thornton 5 
Department of Oceanography 

Naval Postgraduate School 

Monterey, California 93940 


Dr. Taivo Laevastu 2 
EuuToumenval Prediction Agency 

Fleet Numerical Weather Center 

Monterey, California 93940 


Lt. Jose Saldanha 1 
Mastituto Hidrografico 

Ministerio de Marinha 

DisbDoa, Portugal 


A Herrmann 2 
Pevvarics Branch 

U. S. Army Waterways Experiment Station 

Moros os Engineers, PO Box 631 

waeksbure, Mississippi 39160 


Ur. Edward. koehm al 
Ensjneering Division 

5-5! Engineers 

500 N. Los Angeles 

Los Angeles, Calif. 





hee George М, Fisackerly di 
Harbor Entrance Section 

U.S. Army Corps of Engineers 

meer Box 031 

Vicksburg, Mississippi 39180 


Mr. Stuart L. Kufferman 1 
College of Marine Studies 

University of Delaware 

Newark, Delaware 19711 


i. Dennis Fa POLIO E 
College of Marine Studies 

University of Delaware 

Newark, Delaware 19711 


ia. sneldon Μ. Lbazanoff j 
Fleet Numerical Weather Central 
Monterey, California 93910 


E J. A, Galt 1 
Шот отеп of Oceanography ` 

Naval Postgraduate School 

Monterey, California 93910 


EDD. J. Baumgartner I 
Mederal Water Quality Administration 

Borchwest Region 

Pacific Northwest Water Laboratory 

200 Southwest 35th street 

Corvallis, Oregon 9/330 


PER. C. Dean Ш 
Coastal and Oceanographic Engineering Department 
Umiversity ol Florida 

Gainesville, Florida 32601 


environmental and Fishery Forecasting Center T 
National Marine Fisheries Service 

National Oceanie and Atmospheric Administration 
US. Department of Commerce 

C/O Fleet Numerical Weather Central 

Monterey, California 93940 : 


ШИ Ооп НОП 1 
Special Assistant to Director of Navy Laboratory 
tede 034 

Naval Ship System Command 

Кс, D. C. 20350 


136 





EU. 


M. 


e2., 


ES. 


24, 


25. 


26. 


ET. 


BS. 


ES. 


L i Charles A. Farrel 

Ocean services 

Naval Undersea Research and Development Center 
sen Diego, California 92132 


CDR Richard D. Fasing 

Rat ution Control Oftice 
Mica Naval District 

san Diego, California 92132 


Pp Josepn D'Emidio 

OP H5 

I!wreobLropr of ihe Environmental Protection Div. 
Pentagon Building, Room B516 

Mesaington, 0. С. 20350 


ο. απ, b. Reisbig 

Mechanical Engineering School 
Mi versity of Missouri 

Rolla, Mo. 65H01 


Mr. H. б. Ramming 
Institut für Meereskunde 
Universitat Hamburn 
Hamburg, Germany 


Dr. W. Hansen 

Institut fur Meereskunde 
Universitat Hamburg 
Hamburg, Germany 


rector Instátuto Nacional. de Canalizaciones 
e. Andrés Bello 

buco Atysntio 69 

Los Palos Grandes 

Maracas, Venezuela 


раси de Ingenieria 

imsituro Nacional de Canalizaciones 
Av. EL Milagro 

Maracaibo, Venezuela 


Departamento de Oceanografía 
Universidad Simon Bolivar 
tracas, Venezuela 


Departamento de Oceanografia 


Mauversgydad de Oriente 
Cumana, Venezuela 


ү 





B. 


EL. 


ED. 


ES. 


34, 


OR 


3T. 


38. 


29. 


ОЛСО Ор 

National Oceanic and Atmospheric Administration 
ШЕ Department 57. Commerce 

asno DOS 20235 


Division of Oceanography 
Environmental Science Service Administration 
Silver Springs, Marylana 20910 


Commanding Officer 

Office of Naval Research Branch Office 
Box 39 

Fleet Post Office 

New York 09510 


Chief of Naval Research 

Ocean Science and Technology Group 
Code 480 

Office of Naval Research 
Washington, D.C. 20360 


Chief of Naval Research, Code 480 
Peis ον, Ostenso 

Office of Naval Research 
Besshinpton, D.C. 20360 


Oceanographer of the Navy 

Utilice ol the Oceanographer of the Navy 
КОО ОО орг Washington St. 

Alexandria, Virginia 22314 


υπ οσο 

Coastal Engineering Research Center 
temps οἳ Engineers, U.S. Army 

peel bittle Falls Road, N.W. 

Cn econ, D.C. 20315 


Commandant 

We. Coast Guard 
Headquarters 

Ea neton, D.C. 20591 


Erector, Naval Research Laboratory 
meets Lechnical Information Officer 
| Shington, D.C. 20390 


Barector 

Office of Naval Research Branch Office 
1030 East Green St. 

Pasadena, California 91101 


136 





HO. 


41. 


42. 


ae 


HH. 


s Pres Orite 

Comite de Hidrografía j 

Instituto Panamericano de Geografia e Historia 
Ex-Arzobispado 29 

Mexico 18, D.F. Mexico 


Dirección de Hidrografía y Navegación 
Camandancia General de la Marina 
Apartado 6745 

"О оле, Venezuela 


ie Rafael Steer 
лу еш: т Ne 59-15 
ParranguTlla "omnis 


Dr. Noel Boston 

Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 


Com Kari L. Sehriner 


Pear ters, Filth Naval Distrac. 
ШОЛ ОТК Virginia 23511 


2 





PA d 
Secunty Classification 
аР. р т А тте τω” ; 


e 


e OCITA FARAPA Co ALIE N ον PIROT EYE g T OE PILOT AE NL TG IIRL IT NE SPA, ποι, PET LTT TD 
л 


DOCUMENT COHTROL DATA - R & D L 


a 







(Security classification of title, body ol abstract end indexing ennotation must be entered when tie overall report is clasclfled) | 
ra e a oqa σσ A A rs md N 
GINA TING ACTIVITY (Corporate author) 2a, REPORT SECURITY CLASSIFICATION i 


Dnelacesıı jed 


2b. GROUP 













Naval Postgraduate School 
Monterey, California 93940 


ORT TITLE 


Comparison of Numerical and Hydraulic Oceanographic ἐπ "οποιο 
models 


SCRIPTIVE NOTES (Type of report end, inclusivo detes) 


ITHOR:SIı (First name, middle initial, laet namo) 


Leopoldo Salas N. 


PORT DATE 7a. TOT AL NO. OF PAGES Ἢ 7b. NO. OF REFS -A 
| September 1972 ΠΠ 12 


ONTRACT OR GRANT NO. 98. ORIGINATOR'S REPORT NUMBER(S) 
PROJECT NO. 


Cb. OTHER REPORT NO(S) (Any other numbers that may be eselgned 
thie report) 


DISTRIBUTION STATEMENT 


ПКО od for publie release; ein Unico a. 


Barat sua A ων Í —ac лр чечан 
UPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 






Naval Postgraduate School 
Monterey, California 93910 








ıpplied to San Diego Bay, and the results compared both with the 

hydraulic model and the real data obtained by field measurements. 
Phis allows one of the few good comparisons between numerical and 
атаба 1с models for the песо бас иа: СОПОЛ ΠΡ. 

The bay was divided into two sections that were run separately 
in order to obtain the desirable spatial resolution. iS division 
required solving the problems of proper ООП п ас шшр Tecos 
niques between both portions. The solution involved the addition of 
an appended pseudo-bay to the first section of the model in order to 
compensate for the correct tidal prism. Tae effeems oí a proposed 
second open entrance in the southern part of the bay Were studied. 
This resulted in an increase of flushing in the southern Dome iO ππ 
Ehe bay but caused the currents in the center of the bay to be small 
which decreased dispersion in the central portion ОПЕ п Day. 1η 
general, both models produced similar and reliable results, but 
there was a considerable reduction of cost and time with the 
numerical model. 


ABSTRACT 
The Hydro-Numerical Prediction Model developed by Hansen is | x 





Du TS (PAGE 1) 


N 0101-807-6811 ES Security Clessification 


2- 51406 





Security Clessification 





A CIA SFOS DERELICT MEI TRO th МАРЕ ч Т САТУМ TE IO AO RIAL PDA US Uu AIS EM RHET 
LINK A |: LINK B LIHK C k 
KEY WORDS A - 
wr | ποτε] wT [токе] жт {| 
A AI rs aa πετπμάκ. Ey SR > 


Numerical Oceanographic Model 


Hydraulic Model 


Dee A rn 


san Diego Bay 


PETI 








)D 1473 taco 


erga a aaa mama namam ar a ge aa ANNAN 
Z/N 0101-807-6821 n Security Classification A-31409 





^6 AUG 74 a 9 
‚9 Јум» 23452 
a «9 JAM [ὸ 22660 
| 26406 





Thesis 141353 


51526 Salas Römer 
Γς.1 Comparison of numeri- 


cal and hydraulic ocean- 
ographic prediction mod- 
els. 





